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The power output from given sources is usually ascertained via an 
energy flux integral over the normal directions to a remote (far field) 
surface; an alternative procedure, which utilizes an integral that specifies 
the direct rate of working by the source on the resultant field, is described 
and illustrated for both point and continuous source distributions. A com- 
parison between the respective procedures is made in the analysis of sound 
radiated from a periodic dipole source whose axis performs a periodic plane 
angular movement about a fixed direction. Thus, adopting a conventional 
approach, Sretenskii (1956) characterizes the rotating dipole in terms of an 
infinite number of stationary ones along a pair of orthogonal directions in 
the plane and, through the far field representation of the latter, arrives 
at a series development for the instantaneous radiated power, whereas the 
local manner of power calculation dispenses with the equivalent infinite 
aggregate of sources and yields a compact analytical result. 
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ON SOURCE RADIATION 


by 

Harold Levine 

Joint Institute for Aeronautics and Acoustics 
Stanford University 


§1. Introduction 

Determinations of the energy radiated by different source types 
(electromagnetic, acoustic, or elastic) are generally based on the far 
field or asymptotic form of the wave function, although there exists 
an independent way, less well exploited, of calculating the radiation from 
a system: this involves a direct consideration of the rate at which the 

source (electrical current, acoustical force, •*•) does work on the field 
and thus supplies the energy loss by radiation. 

Having regard, specifically, for the microscopic (Maxwell-Lorentz) 
equations 


V X b = 


1 8e ^ 1 

c ^ 


V • e = p 


( 1 ) 


V X e = - 


9b 

at 


b = 0 


of an electromagnetic field described by the two vectors e(r,t) (the 
electric field intensity) and b(r,t) (the magnetic field intensity), with 
the given source functions ^(r, t) (the current density) and p(r, t) (the 
charge density), it is readily shown that 


aw 

at 


V • 



( 2 ) 


where 
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W = 


2 

e 



represents the electromagnetic energy density and 


S 


ce X b 


(3) 


(4) 


characterizes the energy flux (Poynting) vector. An integrated version of 
the local energy balance relation (2) , 

- / J • : d? . If + P , (5) 


associates the rate of working by the source current with the combined 
temporal change of the total electromagnetic energy of the system, 


e 


= f 




dr , 


( 6 ) 


and the net amount of energy radiation therefrom per unit time, or power 
output 

P = / n • ce X b dA , (7) 


reckoned on a distant surface with outward (unit) normal vector n . Thus, 
if the (activity) integral be displayed m the component form (5) , both the 
stored and radiated energy measures are independently and simultaneously 
obtainable. 

Evidently, the electric field within the source region alone has a 
relevance to the aforesaid integral and if the spatial extent of this region 
is sufficiently limited so that the charge/ current density distributions 
undergo small change in the time required for light to traverse the region 
at speed c , the effect of retardation is slight; then the source functions 
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with a retarded time argument which enter into expressions for the field 
intensities can be developed in powers of 1/c , viz.. 


1 1 ^ 1 * |2 9 


p(f,t - ' ' ) = p(r,t) - -^Ir-r' | p(r,t) + -^|r-r' p p(r,t) - 


- — o |r-r* I — ^ p(r,t) + . . . . 

6c 9t 

The leading order estimate for the power radiation deduced from the 
integral in (5) on the basis of such a development proves to be 

P = — ^ (9) 

3c dt 

where p(t) = / r p(r, t) dr denotes the electric dipole moment of the source. 
In the particular circumstance of a dipole moment which oscillates harmonically 
at a single frequency, say 


p(t) = Pq cos u)t , 


the rate of emission of energy has the instantaneous magnitude 


T3 2 w 2 2 

P = ^ 1 — Pn cos 0>t 

- 3 47T 0 

3c 


and fluctuates about an average value 


3 Pq • 


The characterization (9) also obtains after integrating the energy flux over 
a large spherical surface (in the far field), whose radius merely determines 
the time of emission of the observed field; this implies, if the surface 
radius r is centered on a harmonically varying electric dipole, that 


n 2 0) 2 2 . r. 

3c 


with the same average value as in (12) . 
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Having regard, in a similar context, for the linearized system 
of equations which link the pressure and velocity disturbances p(r,t) 
and v(r,t) , in a homogeneous medium with density and sound speed c , 

to the action of a prescribed force F(r , t) , viz., 

^ I? + PQ ^ ^ ° 

c 

and (13) 

Po I = « 


it is readily verified that the local energy balance relation takes the form 



V 




S = F 


V ; 


(14) 


here 


T7 1 2 

" ■ 2 V 




designates the energy density and 


S = pv 


(15) 


(16) 


IS the acoustic energy flux vector. The energy radiation rate from a region 

wherein an external force is applied can thus be inferred from the volume 

integral of F • v , after isolating a time derivative contribution thereto. 

— y 

If the motions are Irrotational and initiated by a source term Q(r,t) 
on the right-hand side of the first equation (13) , then 


and 


V = -V0 , 


P = P, 


n 

0 9t 


(17) 


(18) 
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(where the latter representation follows from the second equation (13)) 
and the velocity potential 0(r, t) satisfies an inhomogeneous wave equation 


c 3t Pq 


(19) 


Multiplication in (19) by 90/ 3t and subsequent rearrangement yields the 
balance equation 



V 


S = Q(^,t) 1^ 


( 20 ) 


with energy density and flux measures. 


and 


W = I Pq(V0)^ + 


2PqC 


- (^)^ 
2 '■gt'' 


( 21 ) 


J . -Po If V0 (22) 

that are the appropriate versions of (15), (16) when (17), (18) hold. Thus, 
^ 90 

the product Q(r,t) serves as a measure of the source activity and its 

0 t 

volume integral can be utilized to determine the net energy radiation. 


The actual manner of calculating energy radiation by analysis of 
integrals over a source region, which this paper aims to describe, commences 
in the next section. Point source models are adopted at first, thereby 
allowing comparison with results otherwise found and, in particular, those 
obtained by Sretenskii (1956) for rotating dipoles. A velocity potential 
and the scalar wave equation underlie the initial presentation, whereas the 
local velocity and a vector wave equation take over the corresponding roles 
in the later account of technique. Finally, some consideration is given to 
extended source distributions. 
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§2. Sound Radxation by Simple and Dipole Sources 

On employing the retarded potential solution of the inhomogeneous 
equation (19) , 


0(?,t) 


Q(r',t - 


r-r' 


4irp^ 


dr' 


I r-r 


(23) 


and assuming that the source function Q(r,t) has a localized nature, i.e.. 


Q(?,t) = 6(?)f(t) , 


(24) 


there emerges the spherical wave function 


0(r,t) 


4iTpQ r 


(25) 


descriptive of a simple source, with variable strength f(t) , at the origin. 
When the latter and its time derivative are expanded in powers of 1/c , it 
follows that 


0(r,t) = T f(t) + *f'(t) + O(^). 


4ttp 


2c 


and 


V (t) + 0(4) ■ 


dt 


4irp 


(26) 


(27) 


0 


with dots symbolizing differentiation. Substitution of (27) into the source 
activity integral [cf. (20)] yields 


I = / Q(^,t) 1^ d? = / 6(?)f(t) 1^ dr 


= / 6(?)f(t) 

= A f^(t) . 

dt 2 


1 ht) 

4ttPq r 

1 j. S(r) 
4iTpQ r 


1 

4ttPqC 

-> 

dr 


f(t) + 0(^)- d? 


1 

4-rrpQC 


c 

f(t)f(t) , 


(28) 
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and, inasmuch as 

f(t)f(t) = (f(t)ht)) - (f(t))^ , 

It can be deduced that- 

^ “ <29) 

expresses the positive definite rate of energy radiation. If the source 
strength varies harmonically, say 


f (t) = m cos oot , 


then 


p = ^ 31„2 o.t 
47rppC 


(30) 


with an average value 


m2 2 
p = ^ m 




(31) 


The fact that the time derivative term in (28) contains a divergent integral, 
r 5(r) 

J — - — dr , IS attributable to the idealization of a point source. 

Let the point source function (24) be replaced by another, 


Q(r,t) = 5(r)f (t) , 


(32) 


which befits a dipole oriented in the x-direction; then 
, 13 

2 _ 


4«Po 

3x \ 

r 

1 


f(t) 

4>rPo 

3x 3 

r 

1 

4>iPo 

|f(t) 

9 

9x 


2c 


6c~ 


3^ 

r_ 

^4 

" 7 


<i) +-i 


• • v • 'V't* 

y f(t) - f(t) + 0(^) 

2rc 3c c 


02 ) 
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and 


90(r,t) 

9t 


— 1 — -ff(t) — (— ) 

4ttPq 3x 


■V- •«» V 

+ f(t) - f(t) + 


2rc 


3c' 


0(^)| 

r 


so that 


I = / Q(r, t) -|^ dr = / 6(r)f(t) dr 




h d^l 

^ J 2c^ 


f(t)f(t) / 5(r) (^) dr 


12ttPqC' 


f(t)f(t) . 


After invoking the determination 

2 2 

/ 6(r) (-^) dr = / 6(x)6(y)6(z) ^ ^ dx dy dz 

r 


2 . 2 

lim 

y , z^O ; x-^0 r 


= 0 


and the rearrangement 


f(t)*f(t) =-^ [f(t)f(t) - f(t)f(t)] + (f(t))^ 


the power output of the dipole is identified as 


(34) 


1 •• 2 

P = ^ (f(t))^ . 

12ttPqC 

Thus, if the dipole moment varies harmonically, and f 
pressions 

4 2 . 

0) m 2 ^ 

P = ;r cos (jJt 


and 


12itPqC' 


4 2 

p = tu m 


24itPqC' 


(35) 

m cos cot , the ex- 

(36) 

(37) 


obtain. 
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Consider next a dipole with harmonically varying moment, whose axis 
rotates at a uniform angular velocity in the x,y plane; the affiliated wave 
function is 


with 


and 


4iTpQ 1 3x ( r ) 3y 1 r 


f^(t) = m cos fit cos wt = Y m{cos ((iH-fi) t + cos(to-fi)t} 


f^Ct) = m sin Qt cos cot = m{sin((vff^) t - sin(co-^^) t} 


(38) 


(39) 


where o) , f2 designate the natural and precession frequencies, respectively. 

On utilizing the expansions for f^ , f^ in (38) which feature rising powers 
of r/c and the concomitant source function 


Q(r,t) = {f^(t) -^ + f^Ct) -^}6(r) 


(40) 


It can be verified that 


I = / Q(?,t) If = +P 


with 


^ = si: ^ 7: A (- 

^0 I 9x dy 




12ttPqC' 


+ f2(t)f2(t) - f2(t)f2(t); 


and 


P = {(fj^Ct))^ + (f^Ct))^} . 

12ttPqC 


( 41 ) 
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Since 


f^(t) = ” Y cos(oi>ffi)t + (a>-f2)^ cos((A>-J^) 1 1 

and 

f^Ct) = sin(CiH-f^)t - (co-f^)^ sin(W“^2)tj> 


the power radiated by the rotating dipole has the instantaneous and time 
average magnitudes, 


p = 1 — J cos 2ootl 

12ttPqC^ ^ i ^ 


m 


24irpQc' 


^ (a 3^ + 6oi)^n^ 4- cos 2oot^ 


and 


P . 


(42) 


(43) 


respectively, the latter in agreement with the prediction of Sretenskii, 

If the axis of the dipole oscillates symmetrically, within a limited 
angular sector, about the x-direction, the characteristic source functions 
in (38) are expressed by 


and 


f^(t) = ra cos ti3t cos(Y cos ^^t) 


f^(t) = m cos tJ^t sin(Y cos ^t) 


(44) 


Straightforward differentiation of these functions , as called for in the 
power representation (41) , yields the explicit result 
2 


P = 


m 


127TP^c' 


jcos^ wt(co^+Y^^ sin^ 


2 2 1 
+ (2Y^ sin cot sin fit - Y^ cos cot cos fit) j 
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12ttPqC 


f 0 / 0 O 0 0 0/ o // / 

cos 03t((ja + 2y (a) sin fit + y cos + y sin f2t) 


+ sin^ 0)t sin"^ fit - y^^cofi^ sin 2o3t sin 2fit j 


wherein the individual magnitudes of o) and fi can be arbitrarily set. 

Sretenskii, whose power analysis utilizes both the far field measure 
of the dipole source and the development 

00 

ly cos fit „ n^ . V infit 
e ’ = Z 1 J (y)e (46) 

n=-oo 

(containing integer order Bessel functions J ), arrives at the result 

n 


where 


12ttPqC 


3 [r^(T) + r^Cx)] , ^ t - I 


“ 2 ~2 

r^(x) + ir 2 (t) = E e cos(arfnr2)x 


He defines, furthermore, an average value of P , viz. 


P = Lim (I / P(t) dtl , 

T-K» \ 0 J 


and, subject to the proviso that no relationship of the form 2(jo/fi = N 
(integral) exists, obtains the expression 

2 00 

P = — 3 Z ((Vfnn)‘^j^(y) 

24ttPqC n=-oo ^ 

which, on the basis of Bessel summation, acquires its ultimate version 


24ttPqC 


3 + | Y^)}. 


The power representation (45), which has finitely many terms in contrast 
with (47) - (48), lends itself to integration and subsequent averaging over any 


span of time, say 0 < t < T ; specifically, 
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T 

/P dt 


m 


0 24ttPqC' 


{t(o3^ + + yV(|- + | yh) 


M 1 2 2 1 2 ^ 3 4 ^ 

+ (-j- - 2 Y ^ Y — + Ye Y -7;^) sin 2o)T 


qqo TOQ *3 / *3 

+ (--^ybi^ + -^yil - Yg Y ^ 2f2T 


, ,1 2 2^2 ,1 2 3 , 1 2 4 5 4 4. sin 2(a^ff])T 

+ (7Y0)fi +7Y0)J^ + 8 ^^ 


^ A 2 2^2 12 3^ 1 2 4 5 4 4 sin 2(o)-fi)T 

+ (— YO)fi --^ycoQ +-gY^ 


16 


03-^2 


+ 


/ sin 2(gyf2S7)T sin 2 (03-2^2) T 
32 \ 0) + 2Q 03 - 2fi 


)} • 


(52) 


Thus, the non-trigonometric terms of (52) dominate in the limit T ->- <» and the 
average value P , as defined by (49) , corresponds exactly with (51) . If the 
span 

T = 2i\/Q. 


equals the precessional period of the dipole, the appertaining average power 

2i\/Q, 


A 

2ty 


f P dt = P 


0 


247TPqC 


r 4 . ^ 2 2^2 ^ 2^4,1^ 3 2, 

^|o3 +3yu3fi + yf^(‘ 2 +'gy) 

^ ,1 3 1 2 _2 ^ 1 2 ^ 3 4 4™ 

+ — (-;ra) - ^ y (i)n + T y — + TFY — ) Sin 

2tt 2 2 ' 4 ' 0) 16 0) 


+ 


1 2 2_2 ^ 1 2 _3 ^ 1 2_4 5 4_4 

7 -Ywfi + — y(^ ~ tz y ^ 

H Z o lo 


(jj + 

1 
8 


1 2 2_2 1 2 _3 ^ 1 2^4 5 4_4-, 

y (j} Q -yYwfi +-gY(^ ~T6^^ 


y 

32 


( 

\ w + 


+ 


to - n 

1 


n 4 ttu) 

2it 


2f^ (0 - 2fi 


\ n . 47T(o ) 


( 53 ) 
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has the limiting estimate 

2 

P = — o (yV(1 + 7 Y^) + O(fiV)} , « » 0) (54) 

24ttPqC-^ *• ^ ^ 

where the explicit (leading) terms are multiplied by a factor of 2 relative 
to the like pair inferred from (51) . 

The average power relative to a period based on the intrinsic frequency 
of the dipole, T = 2i\/b) , can be directly ascertained from (51) and its 
estimate readily secured if 03 » . 


§3. Sound Radiation from Point Forces 

To commence the analysis of power radiation associated with a vectorial 
source function, namely an impressed force, suppose that the latter has a 
variable magnitude and a fixed orientation, i.e., 

F(r,t) = i^(?,t) (55) 

where i,j,k are the trio of unit vectors in the x,y,z directions, respectively. 
Then 


V X 


F = j 


8z 


- k 


3y 


(56) 


and the vorticxty relation which is an immediate consequence of the second of 
the basic equations (13) , namely 


Pq ^ ^ ^ ^ ^ ’ 


(57) 


admits the integral 


V X V 


V _oo — oo * 


(58) 
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After eliminating the pressure p from the system (13) , it follows that 

- Pof] 

c 

and thus, relying on the vector identity 


VxVxv=VV*V“ V^v 


(59) 


1 d ^ 

V v-^^v= - 
c 9t 


(60) 


together with (58), an inhomogeneous wave equation for v obtains 

0 ^ “oo -<» J 

.r| 1 If^-i 

I PqC ^0 -<»\ 3y 8z / J 

+ ji. j. /#dt| . 

I Pq -cc / I pQ -CO 9^9^ J 


Only the x-component of this equation needs to be reckoned with in analysing 

-y- ->■ 

the local source activity or scalar product of F and v , viz. 

.2 


2 _ A_ 

^ 2 2 j^x p 2 9t 


0 


c ^0 -oo \ 9y 9z 


(61) 


c 3t 

If the source function *3^(r,t) has a concentrated spatial nature and is 
temporarily written in the form 

^(r,t) = 5(r) ^ 


the subsequent version of (61) , 

„2 1 8^ \ 

V - — X V = - 

c 9t ^ P^c 


0 


^ 6(r)p(t) +-^p(t){-^ 
^0 V 9y 


+ 


® '«(?) . 
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admits the retarded time solution 


v^(r,t) . 


with the development 


1 2 ( 8 ^ ^ 9 ^ 

3_C^ t ^ \ 3v^ 


AttPqC 


V (r,t) = ■ c^p(t) I 

^ 4ttPqC^ I ^ ^ \ 

-f p(t)( 


3^ ^ 1 


3y^ 3z^ ' ^ 


3y dz / c ^ 

After identifying and grouping all contributions to the integral 


I = / 6(r)p(t)v^(r, t) dr 


which are expressible in the form of a time derivative the remainder, viz. 


P = 


3 3 


1—2 
t(p) . 


AttPqC 

can be interpreted as a measure of power radiation; and, if the replacement 
p(t) = f(t) be made therein, the outcome 


P = 


(f(t))' 


(62) 


12ttPqC' 


is suited to the source function 


«F(r,t) = 6(r)f(t) 


(63) 


with a time varying strength f(t) . 
The pair of functions 


f^(t) = f^ cos fit cos 03t 


and 


(64) 


f (t) = f^ sin fit cos 0)t 

y 0 
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specify instantaneous components of a force with proper frequency w that 
rotates in the x,y plane at the angular frequency ^ about its point of 
application. According to the evident counterpart of (62), 

12ttPqC^ y j 

it follows that 

f ^ 

0 r 2 2 2 2 1 

P = — ^^0) + Q - (03 -fi ) cos 2 03ti (66) 

24ttPqC'^ ^ 

and the average radiated power, 

f ^ 

— -^0 2 2 

P = ^ (o)^+fi^) , (67) 

247TPqC'^ 

corresponds precisely with an integral of the far field energy flux (whose 
expression can be found in Morfey and Tanna (1971)). 

A time harmonic point force whose line of action oscillates through 
a plane angular sector is characterized by the components 


and 


f == cos 03 t cos(y cos Qt) 

X u 

( 68 ) 


f 

y 


f^ cos o)t sin(Y cos fit) . 


On substituting their first derivatives into (65) the radiated power is 
found to have the instantaneous magnitude 


P 



12ttPqC 


3 



sin 03 1 + 


2n2 

Y fi 


2 

cos 


03t 


. 2 

sin 



( 69 ) 


and Its net accumulation or integral between t = 0 and t = T has the 


explicit form 
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/ P dt = ( 

n V 


- (— + — j T - -J sin 2ul + -l-j- ( JJ ;; 


'fi ( sin 2(aH-^)T , sin 2((o-^2)T 


u) + 


0) - S2 


Hence 


P = Lim / P dt 

T-x» 1^0 J 


0 f 2 ^ 1 2^2) 

3 |w + Y Y ) 


24itPqC 

specifies the long time average power output, and 


^0 1 ^ . 4ira) 

, _ 3 1 2 4 Stt 16TT03 n 

127TPqC '• 


. 4Trg) / 1 
32ir ^ la>fn 


- + — H 

f+n by-'Q, j j 


specifies the average output during a precessional period T = 2n/Q . 
» 0 ) it follows from (72) that 


127TPqC 


3 1 2 '^^1^2 /’ 


and this is nearly double the amount predicted by (71) . 


§4. Extended Source Distributions 

When the pointwise nature of the source function Q(r, t) , manifest 
by the representation (24), is set aside and a series development based on 
powers of 1/c invoked for the appertaining wave function (23), namely 


0(r,t) = 




2c 3t 


|?-r'|2 3^ 


Q(r',t) + 

6c 3t 


. . dr' , 
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the concomitant versxon of the source activity integral becomes 

I = /Q(?,t) 

, Q(r,t) ^Q(r',t) ^ ^ ^ 

= lir ^ dr' - 7 -^ / Q(r,t) Q(r',t) dr dr' 

AttPq 1-4.1 AttPqC 


+ ^ f \r-P |Q(?,t) Q(?',t) dr dP 


SttPqC 


8 t' 


— K f |r-r' |^Q(r, t) Q(r',t) dr dr' + 0(l/c^) 

247rpQc"‘ 3t 

On taking note of the various relations 

f df . i A ; Q(?. , t) . q . (£ : . xi dj d?’ 

\r-r'\ 2 dt 1^4,, 


(75) 


(76) 


(justified by the syirmietric character of the integrand relative to r and r’) 
2 

/ Q(r,t) Q(r',t) dr dr' 

3t 

- dJ df' 


^ / Q(?,t) Q(r',t) dr dr' - (^ / Q(?,t)d?)^ , 


and 


(77) 


3 

/ |r-r' |Q(r,t) — ^ Q(r',t) dr dr' 


./ |f-?'l4(Q(?.'=) 3Q(?4 d. , 3 , q(;..t)^|d?d^' 




- / |r-r' I ^ Q(r',t) dr dr' , 

3t^ 

or (on the basis of the aforementioned symmetry) 
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/ d?d?' 


3f 




(78) 


and, finally, 


r “^1 1 2^ ^ 9 Q(r*,t) 

/ Ir-r'l Q(r,t) — ^ dr dr' 


8t 


= -^ / |?-?'|2|q(?,t) d? d?' - d? dr' 

V 9t 9t 


+ f |-_-. 1 2 9 Q(?,t) 8 q(r'.t) ^ 


(79) 


3t 


it follows that 


3t 


I =4£+p 

dt 


where 
P = 


^ (A / Q(?.„ d?)^ ^ / !?-?■ 1 2 3-9(?,t) 3^Q(?;,t) 

^ttPqC ydt / 24 ttp„c-^ 3t 3t 


+ 0(l/c ) 


and 


8tp„ |;_J.| SliPpCdti ; \^2l 


(80) 


(81) 


The prior representations (29) and (35) , appropriate to the case of a 
point monopole or dipole, are directly recovered from the first and second 
terms of (80), respectively, on deploying the source functions (24) and (32) 
If a distributed though unidirectional force F(r,t) = i^(r,t) acts 
the power characterization analogous to (80) has the form 


P - 


4ttPqC 


f 




J_ / 9V(r’,t) ^ 

120TTp^c^ 3t^ 3t"^ 


( 82 ) 
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